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Abstract —Unlike many complex networks studied in the lit¬ 
erature, social networks rarely exhibit unanimous behavior, or 
consensus. This requires a development of mathematical models 
that are sufficiently simple to be examined and capture, at 
the same time, the complex behavior of real social groups, 
where opinions and actions related to them may form clusters 
of different size. One such model, proposed by Frledkln and 
Johnsen, extends the idea of conventional consensus algorithm 
(also referred to as the iterative opinion pooling) to take into 
account the actors’ prejudices, caused by some exogenous factors 
and leading to disagreement in the final opinions. 

In this paper, we offer a novel multidimensional extension, 
describing the evolution of the agents’ opinions on several topics. 
Unlike the existing models, these topics are interdependent, and 
hence the opinions being formed on these topics are also mutually 
dependent. We rigorous examine stability properties of the pro¬ 
posed model, in particular, convergence of the agents’ opinions. 
Although our model assumes synchronous communication among 
the agents, we show that the same final opinions may be reached 
“on average’’ via asynchronous gossip-based protocols. 


I. Introduction 

A social network is an important and attractive case study 
in the theory of complex networks and multi-agent systems. 
Unlike many natural and man-made complex networks, whose 
cooperative behavior is motivated by the attainment of some 
global coordination among the agents, e.g. consensus, opinions 
of social actors usually disagree and may form irregular 
factions (clusters). We use the term “opinion” to broadly refer 
to individuals’ displayed cognitive orientations to objects (e.g., 
topics or issues); the term includes displayed attitudes (signed 
orientations) and beliefs (subjective certainties). A challenging 
problem is to develop a model of opinion dynamics, admitting 
mathematically rigorous analysis, and yet sufficiently instruc¬ 
tive to capture the main properties of real social networks. 
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The backbone of many mathematical models, explaining the 
clustering of continuous opinions, is the idea of homophily or 
biased assimilation [1]: a social actor readily adopts opinions 
of like-minded individuals (under the assumption that its small 
differences of opinion with others are not evaluated as im¬ 
portant), accepting the more deviant opinions with discretion. 
This principle is prominently manifested by various bounded 
confidence models, where the agents completely ignore the 
opinions outside their confidence intervals [2]-[5]. These 
models demonstrate clustering of opinions, however, their 
rigorous mathematical analysis remains a non-trivial problem; 
it is very difficult, for instance, to predict the structure of 
opinion clusters for a given initial condition. Another possible 
explanation of opinion disagreement is antagonism among 
some pairs of agents, naturally described by negative ties [6]. 
Special dynamics of this type, leading to opinion polarization, 
were addressed in [7]-[ll]. It should be noticed, however, 
that experimental evidence securing the postulate of ubiquitous 
negative interpersonal influences (also known as boomerang 
effects) seems to be currently unavailable. Since the first 
definition of boomerang effects [12], the empirical literature 
has concentrated on the special conditions under which these 
effects might arise. This literature provides no assertion that 
boomerang effects, sometimes observed in dyad systems, are 
non-ignorable in multi-agent networks of social influence. 

It is known that even a network with positive and linear cou¬ 
plings may exhibit persistent disagreement and clustering, if 
its nodes are heterogeneous, e.g. some agents are “informed”, 
that is, influenced by some external signals [13], [14]. One 
of the first models of opinion dynamics, employing such 
a heterogeneity, was suggested in [15]-[17] and henceforth 
is referred to as the Friedkin-Johnsen (FJ) model. The FJ 
model promotes and extends the DeGroot iterative pooling 
scheme [18], taking its origins in French’s “theory of social 
power” [19], [20]. Unlike the DeGroot scheme, where each 
actor updates its opinion based on its own and neighbors’ 
opinions, in the FJ model actors can also factor their initial 
opinions, or prejudices, into every iteration of opinion. In 
other words, some of the agents are stubborn in the sense 
that they never forget their prejudices, and thus remain persis¬ 
tently influenced by exogenous conditions under which those 
prejudices were formed [15], [16]. In the recent papers [21], 
[22] a sufficient condition for stability of the FJ model was 
obtained. Furthermore, although the original FJ model is based 
on synchronous communication, in [21], [22] its “lazy” version 
was proposed. This version is based on asynchronous gossip 
influence and provides the same steady opinion on average, no 
matter if one considers the probabilistic average (that is, the 
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expectation) or time-average (the solution Ces^o mean). The 
FJ model and its gossip modification are intimately related to 
the PageRank computation algorithms [22]-[28]. In special 
cases, the FJ model has been given a game-theoretic [29] 
and an electric interpretation [30]. Similar dynamics arise in 
Leontief economic models [31] and some protocols for multi¬ 
agent coordination [32]. Further extensions of the FJ model 
are discussed in the recent papers [33], [34]. 

Whereas many of the aforementioned models of opinion 
dynamics deal with scalar opinions, we deal with influence 
that may modify opinions on several topics, which makes it 
natural to consider vector-valued opinions [4], [35]-[37]; each 
opinion vector in such a model is constituted by d > 1 topic- 
specific scalar opinions. A corresponding multidimensional 
extension has been also suggested for the FJ model [17], [34]. 
However, these extensions assumed that opinions’ dimensions 
are independent, that is, agents’ attitudes to each specific topic 
evolve as if the other dimensions did not exist. In contrast, if 
each opinion vector is constituted by an agent’s opinions on 
several interdependent issues, then the dynamics of the topic- 
specific opinions are entangled. It has long been recognized 
that such interdependence may exist and is important. A set 
of interdependent positions on multiple issues is referred to 
as schema in psychology, ideology in political science, and 
culture in sociology and social anthropology; scientists more 
often use the terms paradigm and doctrine. Converse in his 
seminal paper [38] defined a belief system as a “configuration 
of ideas and attitudes in which elements are bound together 
by some form of constraints of functional interdependence”. 
All these closely related concepts share the common idea of 
an interdependent set of cognitive orientations to objects. 

The main contribution of this paper is a novel multidi¬ 
mensional extension of the FJ model, which describes the 
dynamics of vector-valued opinions, representing individuals’ 
positions on several interdependent issues. This extension, 
describing the evolution of a belief system, cannot be obtained 
by a replication of the scalar FJ model on each issue. For 
both classical and extended FJ models we obtain necessary 
and sufficient conditions of stability and convergence. We also 
develop a randomized asynchronous protocol, which provides 
convergence to the same steady opinion vector as the original 
deterministic dynamics on average. This paper significantly 
extends results of the paper [39], which deals with a special 
case of the FJ model [15], [21], [22] satisfying the “coupling 
condition”. This condition restricts the agent’s susceptibility 
to neighbors’ opinions to coincide with its self-weight. 

The paper is organized as follows. Section II introduces 
some concepts and notation to be used throughout the paper. 
In Section III we introduce the scalar FJ model and related 
concepts; its stability and convergence properties are studied 
in Section IV. A novel multidimensional model of opinion 
dynamics is presented in Section V. Section VI offers an 
asynchronous randomized model of opinion dynamics, that 
is equivalent to the deterministic model on average. We 
illustrate the results by numerical experiments in Section VII. 
In Section VIII we discuss approaches to the estimation of the 
multi-issues dependencies from experimental data. Proofs are 
collected in Section IX. Section X concludes the paper. 


II. Preliminaries and Notation 
Given two integers m and n > m, let m : n denote the 
set {m,m -f 1,. .. ,n}. Given a finite set V, its cardinality 
is denoted by |V|. We denote matrices with capital letters 
A = (oy), using lower case letters for vectors and scalar 
entries. The symbol denotes the column vector of ones 
(1,1,..., 1)^ S K", and is the identity matrix of size n. 

Given a square matrix A = {aij)2j=i, let diagA = 
diag(aii, 022 , • ■ •, Onn) € stand for its main diagonal 

and p{A) be its spectral radius. The matrix A is Schur stable 
if p{A) < 1. The matrix A is row-stochastic if Oy > 0 
and = 1 '^*- Given a pair of matrices A G R™^", 

B G RP^'?, their Kronecker product [40], [41] is defined by 


aiiB 

012.8 • • 

■ 0,inB 

a2iB 

022 ^? 

■ a2nB 


^m2-B * • 

^mnB 


A (directed) graph is a pair Q = (V, £), where V stands for 
the finite set of nodes or vertices and f C V x V is the set of 
arcs or edges. A sequence i = iq ii ir = i' is 

called a walk from i to i'; the node i' is reachable from the 
node i if at least one walk leads from i to f. The graph is 
strongly connected if each node is reachable from any other 
node. Unless otherwise stated, we assume that nodes of each 
graph are indexed from 1 to n = |V|, so that V = 1 : n. 

III. The FJ and DeGroot Models 

Consider a community of n social actors (or agents) indexed 
1 through n, and let x = (xi,..., Xn)^ stand for the column 
vector of their scalar opinions Xi G R. The Friedkin-Johnsen 
(FJ) model of opinions evolution [15]-[17] is determined by 
two matrices, that is a row-stochastic matrix of interpersonal 
influences W G R”^" and a diagonal matrix of actors’ 
susceptibilities to the social influence 0 < A < /„ (we follow 
the notations from [21], [22]). At each stage A: = 0,1,... of 
the influence process the agents’ opinions evolve as follows 

x{k-\-1) = AWx{k)-\-(I — A)u, x(0) = u. (1) 

The values Ui = Xi{0) are referred to as the agents prejudices. 

The model (1) naturally extends DeGroot’s iterative scheme 
of opinion pooling [18] where A = /„. Similar to DeGroot’s 
model, it assumes an averaging (convex combination) mecha¬ 
nism of information integration. Each agent i allocates weights 
to the displayed opinions of others under the constraint of an 
ongoing allocation of weight to the agent’s initial opinion. 
The natural and intensively investigated special case of this 
model assumes the “coupling condition” = 1 — wu Vi 
(that is, A = / — diag W). Under this assumption, the self¬ 
weight Wu plays a special role, considered to be a measure 
of stubborness or closure of the ith agent to interpersonal 
influence. If wu = 1 and thus Wij = OVj f i, then it is 
maximally stubborn and completely ignores opinions of its 
neighbors. Conversely, if wu = 0 (and thus its susceptibility 
is maximal \u = 1), then the agent is completely open to 
interpersonal influence, attaches no weight to its own opinion 
(and thus forgets its initial conditions), relying fully on others’ 
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opinions. The susceptibility of the ith agent = 1 — wu 
varies between 0 and 1, where the extremal values correspond 
respectively to maximally stubborn and open-minded agents. 
From its inception, the usefulness of this special case has been 
empirically assessed with different measures of opinion and 
alternative measurement models of the interpersonal influence 
matrix W [15]-[17], [42]. 

In this section, we consider dynamics of (1) in the general 
case, where the diagonal susceptibility matrix 0 < A < /„ 
may differ from I — diag W. In the case where wu = 1 and 
hence Wij = 0 as z ^ j, one has 0 :^( 1 ) = a;i(0) = Ui and 
then, via induction on k, one easily gets Xi{k) = Ui for any 
fc = 0,1,..., no matter how Xu is chosen. On the other hand, 
if Xu = 0, then Xi{k) = Ui independent of the weights Wij. 
Henceforth we assume, without loss of generality, that for any 
i G 1 : n one either have Xu = 0 and wu = 1 (entailing that 
Xi{k) = Ui) or Xu < 1 and wu < 1- 

It is convenient to associate the matrix W to the graph 
S[W] = {y,£\W]). The set of nodes V = 1 : n of this 
graph is in one-to-one correspondence with the agents and 
the arcs stand for the inter-personal influences (or ties), that is 
{i,j) G £\W] if and only if Wij > 0. A positive self-influence 
weight Wu > 0 corresponds to the self-loop (i,z). We call 
Q = Q[W] the interaction graph of the social network. 

Example 1: Consider a social network of n = 4 actors, ad¬ 
dressed in [15] and having interpersonal influences as follows 


0.220 

0.120 

0.360 

0.300 

0.147 

0.215 

0.344 

0.294 

0 

0 

1 

0 

0.090 

0.178 

0.446 

0.286 


Fig. 1 illustrates the corresponding interaction graph. 



Fig. 1; Interaction graph (/[VF], corresponding to matrix (2) 

In this section we are primarily interested in convergence 
of the FJ model to a stationary point (if such a point exists). 

Definition 1: (Convergence). The FJ model (1) is conver¬ 
gent, if for any vector u G M" the sequence x{k) has a limit 

x' = lim x{k) => x' = AWx' + (f — A)m. (3) 

k—^oc 

It should be noticed that the limit value x' = x'{u) in 
general depends on the initial condition a::(0) = u. A special 
situation where any solution converges to the same equilibrium 
is the exponential stability of the linear system (1), which 
means that AW is a Schur stable matrix: p{AW) < 1. A 


stable FJ model is convergent, and the only stationary point is 

OO 

x' = ^(AIF)'=(/ - A)u ={I- AW)-\I - A)u. (4) 

k=0 

As will be shown, the class of convergent FJ models is in fact 
wider than that of stable ones. This is not surprising since, for 
instance, the classical DeGroot model [18] where A = /„ is 
never stable, yet converges to a consensus value (x^ = ... = 
x'„) whenever W is stochastic indecomposable aperiodic (SIA) 
[43], for instance, is positive for some m > 0 (i.e. W 
is primitive) [18]. In fact, any unstable FJ model contains a 
subgroup of agents whose opinions obey the DeGroot model, 
being independent on the remaining network. To formulate the 
corresponding results, we introduce the following definition. 

Definition 2: (Stubborness and obUvion). We call the ith 
agent stubborn if Xu < 1 and totally stubborn if Xu = 0. 
An agent that is neither stubborn nor influenced by a stubborn 
agent (connected to some stubborn agent by a walk in the 
interaction graph (7[IF]) is called oblivious. 

Example 2: Consider the FJ model (1), where W is from 
(2) and A = I — diag W. It should be noticed that this model 
was reconstructed from real data, obtained in experiments with 
a small group of individuals, following the method proposed 
in [15]. Figure 2 illustrates the graph of the coupling matrix 
AW and the constant “input” (prejudice) u. In this model the 
agent 3 (drawn in red) is totally stubborn, and the three agents 
1, 2 and 4 are stubborn. Hence, there are no oblivious agents in 
this model. As will be shown in the next section (Theorem 1), 
the absence of oblivious agents implies stability. 



U2 


Fig. 2: The structure of couplings among the agents and “in¬ 
puts” for the FJ model with W from (2) and A = I — diag W 

The prejudices m are considered to be formed by some 
exogenous conditions [15], and the agent’s stubborness can 
be considered as their ongoing influence. A totally stubborn 
agent remains affected by those external “cues” and ignores 
the others’ opinions, so its opinion is unchanged Xi{k) = Ui. 
Stubborn agents, being not completely “open-minded”, never 
forget their prejudices and factor them into every iteration of 
opinion. Non-stubborn agents are not “anchored” to their own 
prejudices yet be influenced by the others’ prejudices via com¬ 
munication (such an influence corresponds to a walk in G\W] 
from the agent to some stubborn individual), such individuals 
can be considered as “implicitly stubborn”. Unlike them, for 
an oblivious agent the prejudice does not affect any stage of 
the opinion iteration, except for the first one. The dynamics of 
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oblivious agents thus depend on the “prehistory” of the social 
network only through the initial condition a;(0) = u. 

After renumbering the agents, we assume that stubborn 
agents and agents, influenced by them, are numbered 1 through 
n' <n and oblivious agents (if they exist) have indices from 
n' + 1 to n. For the oblivious agent i we have Xu = 1 and 
Wij = OVj < n'. Indeed, if Wij > 0 for some j < n', then 
the ith agent is connected by a walk to some stubborn agent 
via agent j and hence is not oblivious. The matrices W, A and 
vectors x{k) are therefore decomposed as follows 


W = 


0 


H^12' 
^22 : 


0 



x^(k) 

x\k) 


(5) 


where x^ G R" and and A^^ have dimensions n' x n'. 
If n' = n then x‘^{k), and are absent, otherwise 
the oblivious agents obey the conventional DeGroot dynamics 
x‘^{k + 1) = W‘^‘^x^{k), being independent on the remaining 
agents. If the FJ model is convergent, then the limit = 
lim obviously exists, in other words, the matrix 

k—^oo 

is regular in the sense of [44, Ch.XIII, §7]. 

Definition 3: (Regularity) A matrix A G is called 

regular [44] if a limit A* = lim exists. A regular row- 

k—^oo 

stochastic matrix A is called/uZZy regular [44] or SIA [43] if 
all rows of A* are identical, i.e. A* = IdV^, where v G 
In the literature, the regularity is usually defined for non¬ 
negative matrices [44], but in this paper we use this term for a 
general matrix. In the Appendix we examine some properties 
of stochastic regular matrices, which play an important role 
in the convergence properties of the FJ model. 


IV. Stability and convergence oe the FJ model 

The main contribution of this section is the following 
criterion for the convergence of the FJ model, which employs 
the decomposition (5). 

Theorem 1: (Stability and convergence) The matrix 
is Schur stable. The system (1) is stable if and only 
if there are no oblivious agents, that is, AIF = A^^IF^^. The 
FJ model with oblivious agents is convergent if and only if 
is regular, i.e. the limit = lim exists. In 

k—>oo 

this case, the limiting opinion x' = lim x{k) is given by 

k—¥C!0 



■(/-A“lF“)-i O' 


1 

> 

AlllFl2p^22- 

X = 

0 I 


0 

IFf 


( 6 ) 

An important consequence of Theorem 1 is the stability of 
the FJ model, whose interaction graph is strongly connected 
(or, equivalently, the matrix W is irreducible [44]). 

Corollary 1: If the interaction graph Q[W] is strongly 
connected and A ^ / (i.e. at least one stubborn agent exists), 
then the FJ model (1) is stable. 

Proof: The strong connectivity implies that each agent 
is either stubborn or connected by a walk to any of stubborn 
agents; hence, there are no oblivious agents. ■ 

Theorem 1 also implies that the FJ model is featured by the 
following property. For a general system with constant input 

x{k + 1) = Ax{k) + Bu, (7) 


the regularity of the matrix A is a necessary and sufficient 
condition for convergence if Bu = 0, since x{k) = A^'a;(0) — 
A*a;(0). For Bu f 0, regularity is not sufficient for the 
existence of a limit lim x(k): a trivial counterexample is 

k—foo 

A = I. Iterating the equation (7) with regular A, one obtains 

k oo 

x{k) = A^x{0)+y" A^Bu -^ A*x(0)+y A'^'Su, (8) 

^' fc—>-oo ^' 


where the convergence takes place if and only if the series in 
the right-hand side converge. The convergence criterion from 
Theorem 1 implies that for the FJ model (1) with A = KW 
and B = I — A the regularity of A is necessary and sufficient 
for convergence [34]; for any convergent FJ model (8) holds. 

Corollary 2: The FJ model (1) is convergent if and only if 
A = AW is regular. If this holds, the limit of powers A* is 


A* = lim {AWy = 

k—foo 


0 {I 

0 1^22 


(9) 


and the series from (8) (with B = I — A) converge to 


y(AlF)'=(/-A)u = 

fc=0 


(/- 

0 



Due to (8), the final opinion x' from (6) decomposes into 


x'= A^uA^{AWf{I - A)u. (11) 

fe =0 

Proof: Theorem 1 implies that the matrix A = AW is 
decomposed as follows 

rA“lF“ A11FF12- 
^ - 0 JF22 ’ 


where the submatrix A^^IF^^ is Schur stable. It is obvious that 
A is not regular unless W^^ is regular, since A^ contains the 
right-bottom block (IF^^)^. A straightforward computation 
shows that if IF^^ is regular, then (9) and (10) hold, in 
particular, A is regular as well. ■ 

Note that the first equality in (4) in general/ai/i for unstable 
yet convergent FJ model, even though the series (10) converges 
to a stationary point of the system (1) (the second equality in 
(4) makes no sense as / — AIF is not invertible). Unlike the 
stable case, in the presence of oblivious agents the FJ model 
has multiple stationary points for the same vector of prejudices 
u\ the opinions x{k) and the series (10) converge to distinct 
stationary points unless IF^^u^ = 0. 

As shown in the Appendix, for a regular row-stochastic 
matrix A the limit A* equals to 


A* = lim A^’ = lim (/ — aA) ^(1 — a). (12) 

k—¥OQ OL—fl 

Theorem 1, combined with (12), entails the following im¬ 
portant approximation result. Along with the FJ model (1), 
consider the following “stubborn” approximation 


Xa{k1) = aAWxafk){I — aA)u^ Xa(0) = u, (13) 


where a G (0;1). Hence aA < /, which implies that all 
agents in the model (13) are stubborn, the model (13) is stable. 
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converging to the stationary opinion Xa{k) - > x'^ = {I — 

k—¥oc 

aAW)~^{I — aA)u. It is obvious that Xa{k) -)• x{k) for 

a —>^1 

any fc = 1, 2,..a question arises if such a convergence takes 
place for k = oo, that is, x' as a —>■ 1. A straightforward 

computation, using (12) for A = and (6), shows that this 
is the case whenever the original model (1) is convergent. 
Moreover, the convergence is uniform in u, provided that u 
varies in some compact set. In this sense any convergent FJ 
model can be approximated with the models, where all of 
the agents are stubborn (A < /„). The proof of (12) in the 
Appendix allows to get explicit estimates for \\x'^ — x'\\ that, 
however, do not appear useful for the subsequent analysis. 

V. A MULTIDIMENSIONAL EXTENSION OE THE FJ MODEL 

In this section, we propose an extension of the FJ model, 
dealing with vector opinions xi{k),... ,x„{k) € M*”. The 
elements of each vector Xi{k) = {xj(k),..., x™(k)) stand 
for the opinions of the Ah agent on m different issues. 

A. Opinions on independent issues 

In the simplest situation where agents communicate on m 
completely unrelated issues, it is natural to assume that the 
particular issues x{(fe ),... ,xi^{k) satisfy the FJ model 
(1) for any j = 1,..., m, and therefore 

n 

Xi{k+1) = Xii'^WijXj{k) + {l-Xii)ui, Ui :=Xi(0). (14) 
J=i 

Example 3: Consider the FJ model (14) with W from (2) 
and A = I — diag W. Unlike Example 2, now the opinions 
Xj(k) are two-dimensional, that is, m = 2 and Xj{k) = 
{xj{k),x'j{k))^ represent the opinions on two independent 
topics (a) and (b). The structure of the system, consisting 
of two copies of the usual FJ model (1), is illustrated by 
Figure 3. Since the topic-specific opinions xj(fc), x^(fc) evolve 
independently, their limits can be calculated independently, 
applying (4) to m* = (x*i(0),x^(0),x^(0 ),x|(0))T, i = 1,2. 
For instance, choosing the initial condition 

x(0) = it=[25^ 25,^15, 75, -50, 8^ (15) 

U 2 — X 2 { 0 ) 143=313(0) U 4 =X 4 ( 0 ) 

the vector of steady agents’ opinion is 

x' = [60, -19.3,60, -21.5, 75, -50, 75, -23.2]^. (16) 

B. Interdependent issues: a belief system’s dynamics 

Dealing with opinions on interdependent topics, the opin¬ 
ions being formed on one topic are influenced by the opinions 
held on some of the other topics, so that the topic-specific 
opinions are entangled. Consider a group of people discussing 
two related topics, e.g. fish (as a part of diet) in general and 
salmon. Salmon is nested in fish. A person disliking fish also 
dislikes salmon. If the influence process changes individuals’ 
attitudes toward fish, say promoting fish as a healthy part of 
a diet, then the door is opened for influences on salmon as a 
part of this diet. If, on the other hand, the influence process 


U3a Ll3b 



Fig. 3: The structure of the two-dimensional model (14) 

changes individuals’ attitudes against fish, say warning that 
fish are now contaminated by toxic chemicals, then the door 
is closed for influences on salmon as part of this diet. 

Adjusting his/her position on one of the interdependent 
issues, an individual might have to adjust the positions on 
several related issues simultaneously in order to maintain the 
belief system’s consistency. Contradictions and other inconsis¬ 
tencies between beliefs, attitudes and ideas trigger tensions and 
discomfort (“cognitive dissonance”) that can be resolved by 
a within-individual (introspective) process. This introspective 
process, studied in cognitive dissonance and cognitive con¬ 
sistency theory, is thought to be an automatic process of the 
human brain, enabling an individual to develop a “coherent” 
system of attitudes and beliefs [45], [46]. 

To the best of the authors’ knowledge, no model describing 
how networks of interpersonal influences may generate belief 
systems is available in the literature. In this section, we make 
the first step towards Ailing this gap and propose a model, 
based on the classical FJ model, that takes issues interde¬ 
pendencies into account. We modify the multidimensional FJ 
model (14) (with Xj{k) G M™) as follows 

n 

Xi {^k 1) — XiiC ^ Xj (A;) “1“ (1 Xii'jtli . (17) 

i=i 

The model (17) inherits the structure of the usual FJ dynamics, 
including the matrix of social influences W and the matrix of 
agents’ susceptibilities A. On each stage of opinion iteration 
the agent i calculates an “average” opinion, being the weighted 
sum 'Yhj'^ij^jik) of its own and its neighbors’ opinions; 
along with the agent’s prejudice Ui it determines the updated 
opinion Xi{k -\- 1). The crucial difference with the FJ model 
is the presence of additional introspective transformation, ad¬ 
justing and mixing the averaged topic-specific opinions. This 
transformation is described by a constant “coupling matrix” 
C G henceforth called the matrix of multi-issues 

dependence structure (MiDS). In the case C = Im the model 
(17) reduces to the usual FJ model (14). 

To clarify the role of the MiDS matrix C, consider for 
the moment a network with star-shape topology where all 
the agents follow a totally stubborn leader, i.e. there exists 
j G {1,2,..., n} such that Xjj = 0 and Wij = 1 = Xu for 
any i j, so that Xi{k -|- 1) = Cuj. The opinion changes in 
this system are movements of the opinions of the followers 
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toward the initial opinions of the leader, and these movements 
are strictly based on the direct influences of the leader. The 
entries of the MiDS matrix govern the relative contributions 
of the leader’s issue-specific opinions to the formation of the 
followers’ opinions. Since x^{k + 1) = then Cpq 

is a contribution of the gth issue of the leader’s opinion to the 
pth issue of the follower’s one. In general, instead of a simple 
leader-follower network we have a group of agents, commu¬ 
nicating on m different issues in accordance with the matrix 
of interpersonal influences W. During such communications, 
the ith agent calculates the average '^jWijXj{k) of its own 
opinion and those displayed by the neighbors. The weight Cpq 
measures the effect of the gth issue of this averaged opinion 
to the pth issue of the updated opinion Xi{k + 1). 

Notice that the origins and roles of matrices W and C in 
the multidimensional model (17) are very different. The matrix 
W is a property of the social network, describing its topology 
and social influence structure, which is henceforth assumed 
to be known (the measurement models for the structural 
matrices K,W are discussed in [15]-[17]). At the same time, 
C expresses the interrelations between different topics of 
interest. It seems reasonable that the MiDS matrix should 
be independent of the social network itself, depending on 
introspective processes, forming an individual’s belief system. 

We proceed with examples, which show that introducing 
the MiDS matrix C can substantially change the opinion 
dynamics. These examples deal with the social network of 
n = 4 actors from [15], having the influence matrix (2) and 
the susceptibility matrix A = / — diag W. Unlike Example 3, 
the agents discuss interdependent topics. 

Example 4: Let the agents discuss two topics, (a) and (b), 
say the attitudes towards fish (as a part of diet) in general and 
salmon. We start from the initial condition (15), which means 
that agents 1 and 2 have modest positive liking for fish and 
salmon; the third (totally stubborn) agent has a strong liking 
for fish, but dislikes salmon; the agent 4 has a strong liking 
for fish and a weak positive liking for salmon. Neglecting 
the issues interdependence (C = I 2 ), the final opinion was 
calculated in Example 3 and is given in (16). 

We now introduce a MiDS matrix, taking into account the 
dependencies between the topics 


Cl 


0.8 0.2 

0.3 0.7 


(18) 


As will be shown in Theorem 2, the opinions converge to 

x'c^ = [39.2,12,39,10.1,75, -50,56,5.3]^. (19) 


Hence, introducing the MiDS matrix Ci from (18), with 
its dominant main diagonal, imposes a substantial drag in 
opinions of the “open-minded” agents 1 and 2. In both cases 
their attitudes toward fish become more positive and those 
toward salmon become less positive, compared to the initial 
values (15). However, in the case of dependent issues their 
attitudes toward salmon do not become negative as they did 
in the case of independence. As for the agent 4, its attitude 
towards salmon under the MiDS matrix (18) becomes even 
more positive, compared to the initial value (15), whereas for 
C = I 2 this attitude becomes strongly negative. 


The difference in behavior of the systems (14) and (17) is 
caused by the presence of additional ties (couplings) between 
the topic-specific opinions, imposed by the MiDS matrix C, 
drawn in Eig. 4 in green. Eor simplicity, we show only three of 
these extra ties; there are also ties between the topic-specific 
opinions lb and 2a, 3a, 4a; 2a and lb, 3b, 4b etc. 


U3a U3b 



Eig. 4: The structure of the two-dimensional EJ model (17) 
with C from (18): emergent extra couplings between topic- 
specific opinions 


In Example 4 the additional ties are positive, bringing the 
topic-specific opinions closer to each other. The requirement 
of consistency of a belief systems may also imply the negative 
couplings between different topics if the multidimensional 
opinion contains attitudes to a pair of contrary issues. 

Example 5: Consider Example 4, replacing the stochastic 
MiDS matrix (18) with the non-positive matrix 


C2 


0.8 - 0.2 
-0.3 0.7 


( 20 ) 


The topics (a) and (b), discussed by the agents, are interrelated 
but opposite, e.g. the agents discuss their attitudes to vege¬ 
tarian and all-meat diets. As will be shown, the system (17) 
remains stable. Starting from the same initial opinion (15), the 
agents’ opinions converge to the final value 


x'c^ = [52.3, -30.9,52.1, -33.3,75, -50,68.4, -33.2]”^. 


Similar to the case of uncoupled topics (Example 3), the agents 
converge on positive attitude to vegetarian diets and negative 
attitude to all-meat diet, following thus the prejudice of the 
totally stubborn individual 3. However, the final opinions in 
fact substantially differ. In the case of decoupled topic-specific 
opinions (C = I 2 ) the agents get stronger positive attitudes to 
vegetarian diets than in the case of negative coupling (C = C 2 ) 
while their negative attitudes against all-meat diets are weaker. 

Remark 1: (Restrictions on the MiDS matrix) In many 
applications the topic-specific opinions may vary only in some 
predefined interval. Eor instance, treating them as certainties 
of belief [47] or subjective probabilities [18], [37], the model 
of opinion evolution should imply that the topic-specific 
opinions belong to [0,1]. Similarly, the agents attitudes, i.e. 
signed orientations towards some issues [34], are often scaled 
to the interval [—1,1]. Such a limitation on the opinions 
can make it natural to choose the MiDS matrix C from a 
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special class. For instance, choosing C row-stochastic, the 
model (17) inherits important property of the FJ model [34]: 
if ul = xlifi) G [a,h\\li,j then xlik) G [a,h\\/i,j for any 
k > 0 (here [a, 6] C K is a given interval, i G 1 : n and 
j G 1 : to). This is easily derived from (17) via induction 
on k: if Xi{k) G [a,b]'^ for k = 0,..., s and any i, then 
Cxi{k) G [a, &]"* and hence a;i(s + 1) G [a, &]™. 

In the case of special a and b the assumption of row- 
stochasticity can be further relaxed. For instance, to keep 
the vector opinions Xi{k) in [0,1]™ whenever Ui G [0,1]™, 
the matrix C should be chosen substochastic: Cij > 0 and 
^ 1 for ^riy m. To provide the invariance of the 
hypercube [—1,1]™, the matrix C should satisfy the condition 

m 

IlCjloo = maxY^ \cij \ < 1. (21) 

i ^^ 

7 = 1 

More generally, via induction on k the following invariance 
property can be proved: if 27 is a convex set and Cx G 27 for 
any x G 27, then 27 is an invariant set for the dynamics (17): 
if Ui = Xi(0) G 27 then Xi{k) G 27 for any k. 

In the next subsection, the problems of stability and con¬ 
vergence of the model (17) are addressed. 


C. Convergence of the multidimensional FJ model 

Similar to (15), the stack vectors of opinions x{k) = 
{xi{k)^,..., Xn{k)^)^ and prejudices u = {uj,... ,u^)^ = 
x(0) can be constructed. The dynamics (17) now becomes 

x{k + 1) = [(AIF) (g) C]x{k) + [(/„ - A) 0 Jm], (22) 

which is a convenient representation of (17) in the matrix form. 

We begin with stability analysis of the model (22). In the 
case when C is row-stochastic the stability conditions remain 
the same as for the initial model (1). However, the model (22) 
remains stable for many non-stochastic matrices, including 
those with exponentially unstable eigenvalues. 

Theorem 2: (Stability) The model (22) is stable (i.e. 
KW 0 C is Schur stable) if and only if p[KW)p{C) < 1. 
If this holds, then the vector of ultimate opinions is 

x'q := lim x(k) = (/m„-AlF 0 C')“^[(/„-A) 0 /m]M. (23) 

k—^oo 

If p{C) = 1, the stability of (22) is equivalent to the stability 
of the FJ model (1), i.e. to the absence of oblivious agents. 

Remark 2: Theorem 2, in particular, guarantees stability 
when the original FJ model (1) is stable (p{AW) < 1) 
and C is row-stochastic or, more generally, satisfies (21). 
These conditions, however, are not necessary for stability; the 
system (22) remains stable whenever p{C) < p{KW)~^. 

In the case where some agents are oblivious, the conver¬ 
gence of the model (22) is not possible unless C is regular, 
that is, the limit C, = lim exists (in particular, p{C) < 

k—foo 

1). As in Theorem 1, we assume that oblivious agents are 
indexed n' -\-l through n and consider the decomposition (5). 

Theorem 3: (Convergence) Let n' < n. The model (22) 
is convergent if and only if C is regular and either C* = 0 or 


Xc = 


P = 


Pu, 




(24) 


is regular. If this holds then x{k) - x'q, where 

^—>■00 

■(/-AiiJFii0C')-i 

0 

0 JFf 0 C, 

By definition, in the case where C* = 0 but the limit 
limfc^oo()F^^)^ does not exist we put VF* = 0 . 

Remark 3: (Extensions) In the model (22) we do not 
assume the interdependencies between the initial topic-specific 
opinions; one may also consider a more general case when 
Xi(0) = Dui and hence x(0) = [/„ 0 D]u, where D is 
a constant m x m matrix. This affects neither stability nor 
convergence conditions, and formulas (23), (24) for x'q remain 
valid, replacing P in the latter equation with 

■(/-Aii)0/™ (A“1F12 if22)0C'C'*77' 

0 JF22 0 C'h.TJ 


P = 


VI. Opinion Dynamics under Gossip-Based 
Communication 

A considerable restriction of the model (22), inherited 
from the original Friedkin-Johnsen model, is the simultaneous 
communication between the agents. That is, on each step the 
actors simultaneously communicate to all of their neighbors. 
This type of communication can hardly be implemented in a 
large-scale social network, since, as was mentioned in [15], 
“it is obvious that interpersonal influences do not occur in 
the simultaneous way and there are complex sequences of 
interpersonal influences in a group". A more realistic opinion 
dynamics can be based on asynchronous gossip-based [48], 
[49] communication, assuming that only two agents interact 
during each step. An asynchronous version of the FJ model 
( 1 ) was proposed in [ 21 ], [ 22 ]. 

The idea of the model from [21], [22] is as follows. On each 
step an arc is randomly sampled with the uniform distribution 
from the interaction graph G[W] = {V,£). If this arc is {i,j), 
then the ith agent updates its opinion in accordance with 

Xi{k-\-l) = hi ((1 - yij)xi{k) -f 'YijXj{k))-\-{l-hi)ui. (25) 

Hence, the new opinion of the agent is a weighted average of 
his/her previous opinion, the prejudice and the neighbor’s pre¬ 
vious opinion. The opinions of other agents remain unchanged 

xi{k-\-1) = xi{k) VI fi. (26) 

The coefficient hi G [0,1] is a measure of the agent 
“obstinacy”. If an arc {i,i) is sampled, then 

Xi{k-\-1) = hiXt{k)-\-{1 - hi)ui. (27) 

The smaller is hi, the more stubborn is the agent, for hi = 0 
it becomes totally stubborn. Conversely, for hi = 1 the agent 
is “open-minded” and forgets its prejudice. The coefficient 
G [0,1] expresses how strong is the influence of the jth 
agent on the /th one. Since the arc (i,j) exists if and only if 
Wij > 0 , one may assume that 7 ^ = 0 whenever Wij = 0 . 

It was shown in [21], [22] that, for stable FJ model with 
A = I — diag W, under proper choice of the coefficients hi 










and jij, the expectation Ex{k) converges to the same steady 
value x' as the Friedkin-Johnsen model and, moreover, the 
process is ergodic in both mean-square and almost sure sense. 
In other words, both probabilistic averages (expectations) and 
time averages (referred to as the Cesdro or Polyak averages) 
of the random opinions converge to the final opinion in the 
FJ model. It should be noticed that opinions themselves are 
not convergent (see numerical simulations below) but oscillate 
around their expected values. In this section, we extend the 
gossip algorithm from [21], [22] to the case where A ^ I — 
diag W and the opinions are multidimensional. 

Let G[W] = {V,£) be the interaction graph of the net¬ 
work. Given two matrices r^,r^ such that > 0 and 

lij + lij < 1 , we consider the following multidimensional 
extension of the algorithm (25), (26). On each step an arc is 
uniformly sampled in the set £. If this arc is (i, j), then agent 
i communicates to agent j and updates its opinion as follows 

Xi{k -p 1 ) = (1 - 7 L - yfj)xi{k) + jljCxj{k) + 72 u*. (28) 

Hence during each interaction the agent’s opinion is aver¬ 
aged with its own prejudice and modified neighbors’ opinion 
Cxj{k). The other opinions remain unchanged as in (26). 

The following theorem shows that under the assumption of 
the stability of the original FJ model (22) and proper choice of 
r^,r^ the model (28), (26) inherits the asymptotic properties 
of the deterministic model ( 22 ). 

Theorem 4: (Ergodicity) Assume that p{kW) < 1, i.e. 
there are no oblivious agents, and C is row-stochastic. Let 
F^ = AW and F^ = (/ — A)W. Then, the limit x* = 
lim Ex(fc) exists and equals to the final opinion (23) of the FJ 

k—¥oo 

model (22), i.e. x* = x'q. The random process x(fc) is almost 
sure ergodic, which means that x{k) —>■ x* with probability 
1 , and LP -ergodic so that E||x(/c) — x*||p- > 0, where 

k—¥oo 

1 ^ 

:=-j^—^x{l). (29) 

1=0 

Both equality x* = x'q and ergodicity remain valid, replacing 
= (7 — A)W with any matrix, such that 0 < < 1 — 7 ^ 7 , 

E”=i 7*7 = 1 - A** and 7 ^ = 0 as (i,j) ^ £. 

As a corollary, we obtain the result from [21], [22], stating 
the equivalence on average between the asynchronous opinion 
dynamics (25), (26) and the scalar FJ model (1). 

Corollary 3: Let di be the out-branch degree of the ith 
node, i.e. the cardinality of the set {j : {i,j) G £}. Consider 
the algorithm (25), (26), where Xi G M, {l — hi)di = l — Xa Vi, 
77 G [0,1] and hi'yij = XuWij whenever i ^ j. Then, the 
limit X* = lim Ex(fc) exists and equals to the steady-state 

k—¥oo 

opinion (4) of the FJ model (1): x* = x'. The random process 
x{k) is almost sure and mean-square ergodic. 

Proof: The algorithm (25), (26) can be considered as a 
special case of (28), (26), where C = 1 , 7 A = hi^ij and 7 ^ = 
\ — hi. Since the values yh have no effect on the dynamics 
(28) with C = 1 , one can, changing yk if necessary, assume 
that F^ = AW. The claim now follows from Theorem 4 since 
1 - 7*7 = hi > 7 A and 1% = (1 “ hi)di = 1 - A^*. ■ 


Hence, the gossip algorithm, proposed in [21], [22] is only 
one element of a broad class of protocols (28) (with C = 1), 
satisfying assumptions of Theorem 4. 

Remark 4: (Random opinions) Whereas the Ces^o- 
Polyak averages x(fc) do converge to their average value x*, 
the random opinions x(A:) themselves do not, exhibiting non¬ 
decaying oscillations around x*, see [21] and the numerical 
simulations in Section VII. As implied by [22, Theorem 1], 
x{k) converges in probability to a random vector Xoo, whose 
distribution is the unique invariant distribution of the dynamics 
(28), (26) and is determined by the triple (A, W, C). 

VH. Numerical experiments 

In this section, we give numerical tests, which illustrate the 
convergence of the “synchronous” multidimensional FJ model 
and its “lazy” gossip version. 

We start with the opinion dynamics of a social network 
with n = 4 actors, the matrix of interpersonal influences 
W from (2) and susceptibility matrix A = I — diag W. In 
Fig. 5 we illustrate the dynamics of opinions in Examples 3- 
5: Fig. 5a shows the case of independent issues C = I 2 , 
Fig. 5b illustrates the model (17) with the stochastic matrix C 
from (18), and Fig. 5c demonstrates the dynamics under the 
MiDS matrix (20). As was discussed in Examples 4 and 5, the 
interdependencies between the topics lead to substantial drags 
in the opinions of the agents 1,2 and 4, compared with the 
case of independent topics. 

It is useful to compare the final opinion of the models just 
considered with the DeGroot-like dynamics' (Eig. 6) where the 
initial opinions and matrices C are the same, however, A = /„. 
In the case of independent issues C = I 2 all the opinions are 
attracted by the stubborn agent’s opinion (Eig. 6a) 

lim x(fc) = [75, -50, 75, -50,75, -50, 75, -50]''^. 

k—¥oo 

In the case of positive ties between topics (Eig. 6b) we have 
lim x(fc) = [25, 25, 25, 25, 25, 25, 25, 25]^. 

k—^oo 

In fact, the stubborn agent 3 constantly averages the issues of 
its opinions so that they reach agreement, all other issues are 
also attracted to this consensus value. In the case of negatively 
coupled topics (Eig. 6c) the topic-specific opinions polarize 

lim x{k) = [65, -65,65, -65,65, -65,65, -65]^. 

k—¥oo 

The next simulation (Eig. 7) deals with the randomized 
gossip-based counterparts of the models from Examples (3) 
and (4). The Ces^o-Polyak averages (Eigs. 7a and 7b) x(fc) of 
the opinions under the gossip-based protocol from Theorem 4 
converge to the same limits as the deterministic model (22) 
(blue circles). However, the random opinions x{k) themselves 
do not converge and exhibit oscillations (Eig. 7c). 

*In the DeGroot model [18] the components of the opinion vectors Xi{k) 
are independent. This corresponds to the case where C = Im- One can 
consider a generalized DeGroot’s model as well, which is a special case of (22) 
with A = In but C ^ Im- This implies the issues interdependency, which 
can make all topic-specific opinions (that is, opinions on different issues) 
converge to the same consensus value or polarize, as shown in Eig. 6. 
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(a) Example 3: topics are independent 


(a) DeGroot’s model: independent issues 




Step, k 


(b) Example 4: topics are positively coupled 


(b) Extended DeGroot’s model: MiDS matrix (18) 




(c) Example 5: topics are negatively coupled (c) Extended DeGroot’s model: MiDS matrix (20) 

Fig. 5: Dynamics of opinions in Examples 3-5. Fig. 6: Opinion dynamics in the extended DeGroot model 


The last example deals with the group of n = 51 agents, 
consisting of a totally stubborn “leader” and iV = 10 groups, 
each containing 5 agents (Fig. 8). In each subgroup a “local 
leader” or “representative” exists, who is the only subgroup 
member influenced from outside. The leader of the first 
subgroup is influenced by the totally stubborn agent, and the 
leader of the ith subgroup (i = 2,, N) is influenced by 
that of the {i — l)th subgroup. All other members in each 
subgroup are influenced by the local leader and by each other, 
as shown in Fig. 8. Notice that each agent has a non-zero self¬ 
weight, but we intentionally do not draw self-loops around the 
nodes in order to make the network structure more clear. We 
simulated the dynamics of the network, assuming that the first 


local leader has the self-weight 0.1 (and assigns the weight 
0.9 to the opinion of the totally stubborn agent), and the other 
local leaders have self-weights 0.5 (assigning the weight 0.5 to 
the leaders of predecessing subgroups). All the weights inside 
the subgroups are chosen randomly in a way that W is row- 
stochastic (we do not provide this matrix here due to space 
limitations). We assume that A = I — diag W and choose the 
MiDS matrix as follows 

_ rO.9 0.1' 

^ ~ [o.l 0.9 ■ 

The initial conditions for the totally stubborn agent are 
a;i(0) = [100,-100]^, the other initial conditions are ran¬ 
domly distributed in [—10,10]. The dynamics of opinions in 
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step, k 

(a) Convergence of the CesMO-Polyak averages, the 
MiDS matrix C = I 2 



Step, k 


(b) Convergence of the Cesaro-Polyak averages, the 
MiDS matrix (18) 



(c) The dynamics of opinions with the MiDS ma¬ 
trix (18) show no convergence 

Fig. 7: Gossip-based counterparts of Examples 3 and 4 

the deterministic model and averaged opinions in the gossip 
model are shown respectively in Figs. 9a and 9b. One can see 
that several clusters of opinions emerge, and the gossip-based 
protocol is equivalent to the deterministic model on average, 
in spite of rather slow convergence. 

VIII. Estimation of the MiDS matrix 
Identification of structures and dynamics in social networks, 
using experimental data, is an emerging area, which is actively 
studied by different research communities, working on statis¬ 
tics [50], physics [51] and signal processing [52]. We assume 



that the structure of social influence, that is, the matrices 
W and A are known. A procedure for their experimental 
identification was discussed in [15], [17]. In this section, we 
focus on the estimation of the MiDS matrix C. As discussed 
in [52], identification of social network is an important step 
from network analysis to network control and “sensing” (that 
is, observing or Altering). 

To estimate C, an experiment can be performed where a 
group of individuals with given matrices A and W commu¬ 
nicates on m interdependent issues. The agents are asked to 
record their initial opinions, constituting the vector u = x(0), 
after which they start to communicate. The agents interact in 
pairs and can be separated from each other; the matrix W 
determines the interaction topology of the network, that is, 
which pairs of agent are able to interact. Two natural types of 
methods, allowing to estimate C, can be referred to as “finite- 
horizon” and “infinite horizon” identification procedures. 

A. Finite-horizon identification procedure 

In the experiment of the first kind the agents are asked to 
accomplish T > 1 full rounds of conversations and record 
their opinions Xj(\)^... ^Xj{T) after each of these rounds, 
which can be grouped into stack vectors a:(l),...,a;(T). After 
collection of this data, C can be estimated as the matrix 
best fitting the equations (22) for 0 < k < T. Given 
a;(0) = u, x(l),..., x{T), consider the optimization problem 

T 

minimize llei ll? subject to 

Ej = x{j) — {AW 0 C x{j — 1) + (/ — A) 0 Jm) . 

As was discussed in Remark 1, to provide the model 
“feasibility”, that is, belonging of the opinions to a given set, it 
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(a) 



-40 I 
-60 


(b) 

Fig. 9: Opinion dynamics of n = 51 agents: (a) the determin¬ 
istic model; (b) averaged opinions in the gossip-based model 


can be natural to restrict the MiDS matrix to some known set 
C G C, which may e.g. consist of all row-stochastic matrices 

m 

'^Cij = l Vi, Cij>0 Vi,j, (31) 

i=i 

or all matrices, satisfying (21). In both of these examples, 
the problem (30) remains a convex QP problem, adding the 
constraint C G C. More generally, one may replace (30) with 
the following convex optimization problem 

minimize /(ei, ■ •., St) subject to 
Ej = x{j) — {AW 0 C x{j — 1) + {I — A) ^ I u) , (32) 

C gC. 

Here / is a convex and positive definite function (that is, 
f{ei,... jEt) > 0 and the inequality is strict unless Ei = 
... = Et = 0) and C stands for a closed convex set. In the 
case where f{Ei,...,ET) = EJ=i kjl or /(ei,---,eT) = 
maxj \Ej\ the problem (32) becomes a standard linear pro¬ 
gramming problem. 

Remark 5: In the case where C is a (non-empty) compact 
set, the minimum in (32) exists whenever / is continuous. 
Moreover, if / stands for a strictly convex norm [53, Section 
8.1] on then the minimum exists even for unbounded 
closed convex sets C since the set of all possible vectors 
{(ei, ..., eT)} is closed and convex, and the optimization 
problem (32) boils down to the projection on this set. For 


instance, the minimum in the unconstrained optimization prob¬ 
lem (30) is always achieved. 


B. Infinite-horizon identification procedure 

The experiment of the second kind is applicable only to 
stable models, which inevitably implies a restriction on the 
MiDS matrix C. As in the previous section, we suppose that 
C G C, where C is convex closed set of matrices. We suppose 
also that all elements of C satisfy Theorem 2, that is, p{C) < 
p{AW)~"^. As discussed in Remark 2, the set C may consist 
e.g. of all stochastic matrices (31) or matrices, satisfying (21). 

The agents are not required to trace the history of their opin¬ 
ions, and their interactions are not limited to any prescribed 
number of rounds. Instead, similar to the experiments from 
[15], the agents interact until their opinions stabilize (“agents 
communicate until consensus or deadlock is reached” [15]). 
In this sense, one may assume that the agents know the final 
opinion x'. The matrix C should satisfy the equation 

x'= AW ®Cx'+ {I - A)®Iu, (33) 

obtained as a limit of (22) as k ^ oo. To find a matrix C best 
fitting (33) we solve an optimization problem similar to (30) 

minimize ||e ||2 subject to 
E = x' - {AW ® C x' -G {In - A) ® Im u) , (34) 

C gC. 

In the case where C stands for the polyhedron of matrices 
(e.g. we are confined to row-stochastic MiDS matrices C), the 
problem (34) boils down to a convex QP problem; replacing 
II ■ II 2 with II • 111 or II • ||oo norm, one gets a LP problem. The 
solution in the optimization problem (34) exists for any closed 
and convex set C for the reasons explained in Remark 5. 

Both types of experiments thus lead to convex optimization 
problems. The advantage of the “finite-horizon” experiment is 
its independence of the system convergence. Also, allocating 
some fixed time for each dyadic interaction (and hence for the 
round of interactions), the data collection can be accomplished 
in known time (linearly depending on T). In many applica¬ 
tions, one is primarily interested in the opinion dynamics on 
a finite interval. This approach, however, requires to store 
the whole trajectory of the system, collecting thus a large 
amount of data (growing as nT) and leads to a larger convex 
optimization problem. The loss of data from one of the agents 
in general requires to restart the experiment. On the other hand, 
the “infinite-horizon” experiment is applicable only to stable 
models, and one cannot predict how fast the convergence will 
be. This experiment does not require agents to trace their 
history and thus reduces the size of the optimization problem. 

C. Transformation of the equality constraints 

Both optimization problems (30) and (34) are featured by 
non-standard linear constraints, involving Kronecker products. 
To avoid Kronecker operations, we transform the constraints 
into a standard form Ax = b, where A is a matrix and x, b 
are vectors. To this end, we perform a vectorization operation. 
Given a matrix M, its vectorization vec M is a column vector 
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obtained by stacking the columns of M, one on top of one 
another [40], e.g. vec (2 ?) = [1, 2, 0,1]^. 

Lemma 1: [40] For any three matrices A, B, C such that 
the product ABC is defined, one has 

vec ABC = (C^ <8) .4) vec B. (35) 


The model (17) with this matrix C gives the opinions 
i(l) = [43.12,14.66,42.54,12.37,75, -50,59.90,7.17]^ 
i(2) = [41.93,13.26,41.73,11.35,75, -50,58.37,6.26]^ 
i(3) = [41.30,12.69,41.12,10.79,75, -50,57.90,5.83]^. 


In particular, for A G and S G one obtains 

vec AB = (/„ 8) -4) vec B — {B^ 8) Im) vec A. (36) 

The constraints in (30), (34) can be simplified. Consider first 
the constraint in (34). Let x[ be the final opinion of the ith 
agent and X' = [a;^,..., x'^ be the matrix constituted by 
them, hence x' = vec4f'. Applying (36) for ,4 = C and B = 
X' entails that [/„ 8) C\x' = [(46'^ 8> Im] vecC, thus [KW 8 
Cja;' = [AIL 8 Im\[In 8 C]x' = [XW{X')^ 8 /m]vecC. 
Denoting c = vec C, the constraint in (34) shapes into 

£ + [KWX'^ 8 Im]c = x' - [{In - A) 8 Im\u, (37) 

where the matrix KWX^ 8 Im and the right-hand side are 
known. The constraints in (30) can be rewritten as 

Ej + [KWX{j - lY(i)Im]c = X{j) - [{In-A)®Im]u. (38) 

Here X{j) is the matrix [a;i(j),..., Xn{j)] and hence x{j) = 
vecX{j). 

D. Numerical examples 

To illustrate the identification procedures, we consider two 
numerical examples. 

Example 6: Consider a social network with the matrix W 
from (2), A = / — diagVF and the prejudice vector (15). 
Unlike Example 4, C is unknown and is to be found from the 
“infinite-horizon” experiment. Suppose that agents were asked 
to compute the final opinion, obtaining 

x' = [ 35 , 11 , 35 , 10 , 75 , - 50 , 53 , 5 ]^. 


Remark 6: Examples 6 and 7, demonstrating the estima¬ 
tion procedures, are constructed as follows. We get the model 
(22) with W from (2), A = / — diaglU and C from (18) and 
slightly change its final value (19) (Example 6) and trajectory 
(Example 7). Due to these perturbations, the estimated MiDS 
matrix does not exactly coincide with (18) yet is close to it. 

IX. Proofs 

We start with the proof of Theorem 1, which requires some 
additional techniques. 

Definition 4: (Substochasticity) A non-negative matrix 
A = (oy) is row-substochastic, if < 1 Vi. Given such 

a matrix sized n x n, we call a subset of indices J C 1 : n 
stochastic if the corresponding submatrix {aij)ij^j is row- 
stochastic, i.e. o-ij = 1 Vi G J. 

OEJ 

The Gerschgorin Disk Theorem [44] implies that any sub¬ 
stochastic matrix A has p{A) < 1. Our aim is to identify 
the class of substochastic matrices with p{A) = 1. As it will 
be shown, such matrices are either row-stochastic or contain 
a row-stochastic submatrix, i.e. has a non-empty stochastic 
subset of indices. 

Lemma 2: Any square substochastic matrix A with 
p{A) = 1 admits a non-empty stochastic subset of indices. 
The union of two stochastic subsets is stochastic again, so that 
the maximal stochastic subset J, exists. Making a permutation 
of indices such that J* = (n' -I- 1) : n, where 0 < n' < n, the 
matrix A is decomposed into upper triangular form 

1 -). 09 ) 


Solving the problem (34), one gets the minimal residual 
||£||2 = 0.9322, which gives the MiDS matrix 


0.7562 0.2438 
0.3032 0.6968 


In accordance with (23), this matrix C corresponds to the 
steady opinion 


x'c = [35.316,11.443,35.092,9.483,75, -50,52.386,4.915]^. 


Example 7: Eor A, W and u from the previous example, 
agents were asked to conduct T = 3 full rounds of conversa¬ 
tion (“finite horizon” experiment), obtaining the opinions 

a;(l) = [42.80,14.05,43.59,12.51,75, -50,61.49,7.18]"^ 
x{2) = [41.31,13.37,41.45,11.43,75,-50,55.48,6.45]^ 
x(3) = [41.74,12.30,40.41,10.84,75,-50,58.99,6.02]^. 


We are interested in finding a row-stochastic matrix C best 
fitting the data. Solving the corresponding QP problem (30) 
with additional constraints (31), the MiDS matrix is 


0.8181 0.1819 
0.2983 0.7017 


where A^^ is a Schur stable n' x n'-matrix (p(A^^) < 1) and 

A22 

is row-stochastic. 

Proof: Thanks to the Perron-Erobenius Theorem [41], 
[44], p{A) = 1 is an eigenvalue of A, corresponding to a non¬ 
negative eigenvector v G M" (here n stands for the dimension 
of A). Without loss of generality, assume that max^ Vi = 1. 
Then we either have Vi = 1„ and hence A is row-stochastic (so 
the claim is obvious), or there exists a non-empty set J C 1 : n 
of such indices i that Vi = 1. We are going to show that J is 
stochastic. Since Vi = 1 for i € J'^ = 1 : n\ J, one has 

1 = tiijVj + Uij < IVi G J. 

Since Vj < 1 as j G J°, the equality holds only if = 
0\/i G J,j ^ J and i-®- <7 is a stochastic set. 

This proves the first claim of Lemma 2. 

Given a stochastic subset J, it is obvious that Oy = 0 
when i G J and j ^ J, since otherwise one would have 
^ Oij > 1. This implies that given two stochastic subsets 

jel:n 

Ji, J 2 and choosing i G Ji, one has ^ + 

JGJ1UJ2 iGJi 
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o,ij = 1- The same holds for i € J 2 , which proves 

iGJ2nJf 

stochasticity of the set Ji U J 2 . This proves the second claim 
of Lemma 2 and the existence of the maximal stochastic subset 
J*, which, after a permutation of indices, becomes as follows 
J* = {n' + 1) : n. Recalling that aij = Q'ii € J^:,j & one 

shows that the matrix is decomposed as (39), where is 
row-stochastic. It remains to show that p{A^^) < 1. Assume, 
on the contrary, that = 1. Applying the first claim 

of Lemma 2 to A^^, one proves the existence of another 
stochastic subset J' C 1 : n', which contradicts the maximality 
of J*. This contradiction shows that A^^ is Schur stable. ■ 
Returning to the FJ model (1), it is easily shown now that 
the maximal stochastic subset of indices of the matrix KW 
consists of indices of oblivious agents. 

Lemma 3: Given a FJ model (1) with the matrix A diag¬ 
onal (where 0 < < 1) and the matrix W row-stochastic, 

the maximal stochastic set of indices J* for the matrix KW is 
constituted by the indices of oblivious agents. In other words, 
j G J* if and only if the jth agent is oblivious. 

Proof: Notice, first, that the set J* consists of oblivious 
agents. Indeed, 1 = KiWij < Xu < 1 for any i & J^, 

and hence none of agents from J* is stubborn. Since = 
OVi G J*,j G (see the proof of Lemma 2), the agents 
from J* are also unaffected by stubborn agents, being thus 
oblivious. Consider the set J of all oblivious agents, which, 
as has been just proved, comprises J*: J 3 J*. By definition, 
Xjj = IVj G J. Furthermore, no walk in the graph from J 
to J'^ exists, and hence Wij = 0 as i G J,j ^ J'^, so that 
'^j^jWij = IVz G J. Therefore, indices of oblivious agents 
constitute a stochastic set J, and hence J Q Therefore 
J = J*, which finishes the proof. ■ 

We are now ready to prove Theorem 1. 

Proof of Theorem 1: Applying Lemma 2 to the matrix 
A = KW, we prove that agents can be re-indexed in a way 
that A is decomposed as (39), where A^^ = is Schur 

stable and A^^ is row-stochastic (if A is Schur stable, then 
A = and and A^^ are absent). Lemma 3 shows that 
indices 1 : n' correspond to stubborn agents and agents they 
influence, whereas indices (n' + 1) : n denumerate oblivious 
agents that are, in particular, not stubborn and hence Xjj = 1 
as j > n' so that = W"^^. This proves the first claim of 
Theorem 1, concerning the Schur stability of K^^W^^. 

By noticing that x^{k) = (W^^)'^x^(O), one shows that 
convergence of the FJ model is possible only when W^^ is 
regular, i.e. (W^^)^ —?► W^^ and hence x^(k) —>■ Wf^u^. If 
this holds, one immediately obtains ( 6 ) since 


p{KW)p{C) < 1. In particular, if p{C) = 1 then the system 
(22) is stable if and only if p[KW) <1. ■ 

The proof of Theorem 3 is similar to that of Theorem 1. 
After renumbering the agents, one can assume that obliv¬ 
ious agents are indexed n' + 1 through n and consider 
the corresponding submatrices W^"^ ^W^^ ,W^^, K^^, used in 
Theorem 1. Then the matrix KW 0 C can also be decomposed 


KW 


V 0 W‘^‘^(^C 


(40) 


where the matrices K^^W^^ 0 C has dimensions ran' x 
mn' and m{n — n') x m{n — n') respectively. We con¬ 
sider the corresponding subdivision of the vectors x{k) = 
x‘^{k)^Y ^nd u = 

Proof of Theorem 3: Since the opinion dynamics of 
oblivious agents is given by x^{k -I- 1) = 0 Cx^(k), 

the convergence implies regularity of the matrix W^^ 0 C. 
The regularity of W"^^ 0 C entails that C is regular. Indeed, 
consider a left eigenvector z of at 1 (that is, zIF = z) 
and denote u = z 0 Im- Since u^(IF^^ 0 C)^ = z ® 
has a limit as fc —>■ cx) and z 7 ^ 0 , C is regular and, in 
particular, p{C) < 1. Obviously, the limit C* = limfc_>oo 
is zero if and only if C is Schur stable, i.e. p{C) < 1. If 
p(C) = 1 then A = 1 the only possible eigenvalue on the 
unit circle {A : |A| = 1}. Consider a right eigenvector z 
of C (possibly, complex), corresponding to this eigenvalue. 
Denoting v = In® z, the matrix (IF^^ 0 C)^v = (IF^^)^ 0 v 
has a limit as fc —>• 00 , and thus IF^^ is regular. The 
necessity part is proved. To prove sufficiency, notice that 
x'^{k) 0 as fc —> 00 (where we put IF* = 0 

when C* = 0, and thus (24) follows from the equation 

x^{k + 1) = [A“IF“ 0 C]x^{k) + [A“IF^^ 0 C]x^{k)+ 

+ [/-A“]0/„m\ 


where A^^IF^^ 0 C is Schur stable. ■ 

To proceed with the proof of Theorem 4, we need some 
extra notation. As for the scalar opinion case in [21], [22] the 
gossip-based protocol (28), (26) shapes into 


x{k + 1) = P{k)x{k) + B{k)u, (41) 

where P{k), B(k) are independent identically distributed 
(i.i.d.) random matrices. If arc {i,j) is sampled, then P{k) = 
A(*’J) and B{k) = B^'X)^ where by definition 

= {Imn - {llj + lljYiel 0 + lljtiel 0 C) , 

= 7^ e.gT ^ 


x^{k + 1 ) = K^^W^^x^{k) + K^^W^^x^{k) + (/ - K^^)u^ 

and K^^W^^ is Schur stable. ■ 

The proof of Theorem 2 follows from the well-known 
property of the Kronecker product. 

Lemma 4: [40, Theorem 13.12] The spectrum of the 

matrix A® B consists of all products XiPj, where Ai,..., A„ 
are eigenvalues of A and pi,, prn are those of B. 

Proof of Theorem 2: Lemma 4 entails that p{KW ®C) = 
p{KW)p{C), hence the system (22) is stable if and only if 


Denoting a := \£\ ^ G (0,1] and noticing that EP(fc) = 
EB(fc) =aY,{i,j)(.s , the follow¬ 

ing equalities are easily obtained 

EP(fc) = Imn - a [Imn “ AIF 0 C] 

EB{k) = a{In - A) 0 Im- 

Proof of Theorem 4: As implied by equations (41) and 
(42), the opinion dynamics obeys the equation 

x{k + 1) = P{k)x{k) + v{k) 


(43) 
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where the matrices P{k) and vectors v{k) are i.i.d. and their 
finite first moments are given by the following 

EP(fc) = (1 — a)I + ahW ® C, Er!(/c) = a{In — A) ® Im u, 

where a G (0;1]. Since P(k) are non-negative, Theorem 1 
from [22] is applicable to (43), entailing that the process x(k) 
is almost sure ergodic and Ej:(fc) —>■ x* as fc —c», where 

= [I - AW 0 - A)® Im]u = x'q. 

To prove the L^-ergodicity, notice that x{k) and x{k) remain 
bounded due to Remark 1, and hence E||a;(fc) — x*||^ —!■ 0 
thanks to the Dominated Convergence Theorem [54]. ■ 

Remark 7: (Convergence rate) For the case of p = 2 
(mean-square ergodicity) there is an elegant estimate for the 
convergence rate [21], [27]: E||a;(fc) — x* IP < xKk + 1), 
where x depends on the spectral radius p{AW) and the vector 
of prejudices u. An analogous estimate can be proved for the 
multidimensional gossip algorithm (28), (26). 

Remark 8: (Relaxation of the stochasticity condition) As 
can be seen from the proof. Theorem 4 retains its validity 
for substochastic matrices, since they are non-negative and 
provide boundedness of the solutions. Furthermore, the proof 
of almost sure ergodicity does not rely on the solutions’ 
boundedness and hence is preserved whenever C is non¬ 
negative and p{C)p{AW) < 1. A closer examination of the 
proof of [22, Theorem 1] shows that it can be extended to the 
case where P{k) with negative entries. The non-negativity of 
C can thus be relaxed, however, this relaxation is beyond the 
scope of this paper. 

X. Conclusion 

In this paper, we propose a novel model of opinion dynamics 
in a social network with static topology. Our model is a 
significant extension of the classical Friedkin-Johnsen model 
[15] to the case where agents’ opinions on two or more 
interdependent topics are being influenced. The extension is 
natural if the agent are communicating on several “logically” 
related topics. In the sociological literature, an interdependent 
set of attitudes and beliefs on multiple issues is referred 
to as an ideological or belief system [38]. A specification 
of the interpersonal influence mechanisms and networks that 
contribute to the formation of ideological-belief systems has 
remained an open problem. 

We establish necessary and sufficient conditions for the 
stability of our model and its convergence, which means 
that opinions converge to finite limit values for any initial 
conditions. We also address the problem of identification of 
the multi-issue interdependence structure. Although our model 
requires the agents to communicate synchronously, we show 
that the same final opinions can be reached by use of a 
decentralized and asynchronous gossip-based protocol. 

Several potential topics of future research are concerned 
with experimental validation of our models for large sets 
of data and investigation of their system-theoretic properties 
such as e.g. robustness and controllability. An important open 
problem is the stability of an extension of the model (17), 
where the agents’ MiDS matrices are heterogeneous. A nu¬ 
merical analysis of a system with two different MiDS matrices 


is available in our recent paper [55], further developing the 
theory of logically constrained belief systems formation. 

References 

[1] P. Dandekai; A. Goel, and D. Lee, “Biased assimilation, homophily, and 
the dynamics of polarization,” PNAS, vol. 110, no. 15, pp. 5791—5796, 
2013. 

[2] R. Hegselmann and U. Krause, “Opinion dynamics and bounded con¬ 
fidence models, analysis, and simulation,” Journal of Artifical Societies 
and Social Simulation (JASSS), vol. 5, no. 3, p. 2, 2002. 

[3] G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch, “Mixing beliefs 
among interacting agents,” Advances in Complex Systems, vol. 3, pp. 
87-98, 2000. 

[4] J. Lorenz, “Continuous opinion dynamics under bounded confidence: a 
suiwey,” Int. J. of Modern Phys. C, vol. 18, no. 12, pp. 1819-1838, 2007. 

[5] V. Blondel, J. Hendrickx, and J. Tsitsiklis, “On Ki'ause’s multiagent con¬ 
sensus model with state-dependent connectivity,” IEEE Trans. Autom. 
Control, vol. 54, no. 11, pp. 2586-2597, 2009. 

[6] A. Flache and M. Macy, “Small worlds and cultural polarization,” 
Journal of Math. Sociology, vol. 35, no. 1-3, pp. 146-176, 2011. 

[7] C. Altafini, “Dynamics of opinion forming in structurally balanced social 
networks,” PLoS ONE, vol. 7, no. 6, p. e38135, 2012. 

[8] -, “Consensus problems on networks with antagonistic interactions,” 

IEEE Trans. Autom. Control, vol. 58, no. 4, pp. 935-946, 2013. 

[9] M. Valcher and P. Misra, “On the consensus and bipartite consensus 
in high-order multi-agent dynamical systems with antagonistic interac¬ 
tions,” Systems Control Letters, vol. 66, pp. 94-103, 2014. 

[10] A. Proskurnikov, A. Matveev, and M. Cao, “Opinion dynamics in social 
networks with hostile camps: Consensus vs. polarization,” IEEE Trans. 
Autom. Control, vol. 61, no. 6, pp. 1524—1536, 2016. 

[11] -, “Consensus and polarization in Altafini’s model with bidirectional 

time-varying network topologies,” in Proceedings of IEEE CDC 2014, 
Los Angeles, CA, 2014, pp. 2112-2117. 

[12] C. Hovland, I. Janis, and H. Kelley, Communication and persuasion. 
New Haven: Yale Univ. Press, 1953. 

[13] W. Xia and M. Cao, “Clustering in diffusively coupled networks,” 
Automatica, vol. 47, no. 11, pp. 2395-2405, 2011. 

[14] D. Aeyels and F. D. Smet, “Cluster formation in a time-vaiying multi¬ 
agent system,” Automatica, vol. 47, no. 11, pp. 2481-2487, 2011. 

[15] N. Friedkin and E. Johnsen, “Social influence networks and opinion 
change,” in Advances in Group Processes, 1999, vol. 16, pp. 1-29. 

[16] N. Friedkin, A Structural Theory of Social Influence. New York: 
Cambridge Univ. Press, 1998. 

[17] N. Friedkin and E. Johnsen, Social Influence Network Theory. New 
York: Cambridge Univ. Press, 2011. 

[18] M. DeGroot, “Reaching a consensus,” Journal of the American Statisti¬ 
cal Association, vol. 69, pp. 118-121, 1974. 

[19] J. French Jr., “A formal theory of social power,” The Physchological 
Review, vol. 63, pp. 181—194, 1956. 

[20] F. Harary, “A criterion for unanimity in French’s theory of social power,” 
in Studies in Social Power, Oxford, England, 1959, pp. 168-182. 

[21] P. Frasca, C. Ravazzi, R. Tempo, and H. Ishii, “Gossips and prejudices: 
Ergodic randomized dynamics in social networks,” in Proc. of IFAC 
NecSys 2013 Workshop, Koblenz, Germany, 2013, pp. 212-219. 

[22] C. Ravazzi, P. Frasca, R. Tempo, and H. Ishii, “Ergodic randomized 
algorithms and dynamics over networks,” IEEE Trans. Control of 
Network Syst, vol. 2, no. 1, pp. 78-87, 2015. 

[23] N. Friedkin and E. Johnsen, “Two steps to obfuscation,” Social Net¬ 
works, vol. 39, pp. 12-13, 2014. 

[24] H. Ishii and R. Tempo, “Distributed randomized algorithms for the 
PageRank computation,” IEEE Trans. Autom. Control, vol. 55, no. 9, 
pp. 1987-2002, 2010. 

[25] -, “The PageRank problem, multi-agent consensus and Web aggre¬ 

gation: A systems and control viewpoint,” IEEE Control Syst. Mag., 
vol. 34, no. 3, pp. 34-53, 2014. 

[26] R. Tempo, G. Calafiore, and F. Dabbene, Randomized Algorithms for 
Analysis and Control of Uncertain Systems. Springer-Verlag, 2013. 

[27] P. Frasca, H. Ishii, C. Ravazzi, and R. Tempo, “Distributed randomized 
algorithms for opinion formation, centrality computation and power 
systems estimation: A tutorial overview,” Europ. J. Control, vol. 24, 
no. 7, pp. 2-13, 2015. 

[28] A. Proskurnikov, R. Tempo, and M. Cao, “PageRank and opinion 
dynamics: Missing links and extensions,” in Proc. of IEEE Conf. Norbert 
Wiener in the 2Ist Century, Melbourne, 2016, pp. 12-17. 


15 


[29] D. Bindel, J. Kleinberg, and S. Oren, “How bad is forming your own 
opinion?” in Proc. of IEEE Symp. on Foundations of Computer Science, 
2011, pp. 57-66. 

[30] J. Ghaderi and R. Srikant, “Opinion dynamics in social networks with 
stubborn agents: Equilibrium and convergence rate,” Automatica, vol. 50, 
no. 12, pp. 3209-3215, 2014. 

[31] M. Milanese, R. Tempo, and A. Vicino, “Optimal eiTor predictors for 
economic models,” Int. J. Syst. Set, vol. 19, no. 7, pp. 1189-1200, 1988. 

[32] Y. Cao, “Consensus of multi-agent systems with state constraints: A 
unified view of opinion dynamics and containment control,” in Proc. of 
American Contol Conference (ACC), 2015, pp. 1440-1444. 

[33] P. Jia, A. Mirtabatabaei, N. Friedkin, and F. Bullo, “Opinion dynamics 
and the evolution of social power in influence networks,” SIAM Review, 
vol. 57, no. 3, pp. 367-397, 2015. 

[34] N. Friedkin, “The problem of social control and coordination of complex 
systems in sociology: A look at the community cleavage problem,” IEEE 
Control Syst. Mag., vol. 35, no. 3, pp. 40-51, 2015. 

[35] S. Fortunato, V. Latora, A. Pluchino, and A. Rapisarda, “Vector opinion 
dynamics in a bounded confidence consensus model,” Int. J. of Modern 
Phys. C, vol. 16, pp. 1535-1551, 2005. 

[36] A. Nedic and B. Touri, “Multi-dimensional Hegselmann-Krause dynam¬ 
ics,” in Proc. of IEEE CDC, Maui, Hawaii, USA, 2012, pp. 68-73. 

[37] L. Li, A. Scaglione, A. Swami, and Q. Zhao, “Consensus, polarization 
and clustering of opinions in social networks,” IEEE J. on Selected Areas 
in Communications, vol. 31, no. 6, pp. 1072-1083, 2013. 

[38] P. Converse, “The nature of belief systems in mass publics,” in Ideology 
and Discontent. New York: Free Press, 1964, pp. 206-261. 

[39] S. Parsegov, A. Proskumikov, R. Tempo, and N. Friedkin, “A new model 
of opinion dynamics for social actors with multiple interdependent 
attitudes and prejudices,” in Proc. of IEEE Conference on Decision and 
Control (CDC), 2015, pp. 3475-3480. 

[40] A. Laub, Matrix analysis for scientists and engineers. Philadelphia, 
PA: SIAM, 2005. 

[41] R. Horn and C. Johnson, Topics in Matrix Analysis. New York: 
Cambridge Univ. Press, 1991. 

[42] C. Childress and N. Friedkin, “Cultural reception and production,” Amer. 
Soc. Rev., vol. 77, pp. 45-68, 2012. 

[43] J. Wolfowitz, “Products of indecomposable, aperiodic, stochastic matri¬ 
ces,” Proceedings of Amer. Math. Soc., vol. 15, pp. 733-737, 1963. 

[44] F. Gantmacher, The Theory of Matrices. AMS Chelsea Publishing, 
2000, vol. 2. 

[45] L. Festinger, A Theory of Cognitive Dissonance. Stanford, CA: Stanford 
Univ. Press, 1957. 

[46] B. Gawronski and F. Strack, Eds., Cognitive consistency: A fundamental 
principle in social cognition. New York, NY: Gulford Press, 2012. 

[47] J. Y Halpern, “The relationship between knowledge, belief, and cer¬ 
tainty,” Annals of Mathematics and Artificial Intelligence, vol. 4, no. 3, 
pp. 301-322, 1991. 

[48] S. Boyd, A. Ghosh, B. Prabhakar, and D. Shah, “Randomized gossip 
algorithms,” IEEE Trans. Inform. Theory, vol. 52, no. 6, pp. 2508-2530, 
2006. 

[49] S. Hedetniemi, S. Hedetniemi, and A. Liestman, “A survey of gossiping 
and broadcasting in communication networks,” Networks, vol. 18, no. 4, 
pp. 319-349, 1988. 

[50] T. Snijders, J. Koskinen, and M. Schweinberger, “Maximum likelihood 
estimation for social network dynamics,” The Annals of Applied Statis¬ 
tics, vol. 4, no. 2, pp. 5214-5240, 2010. 

[51] W.-X. Wang, Y.-C. Lai, C. Grebogi, and J. Ye, “Network reconstruction 
based on evolutionaiy-game data via compressive sensing,” Phys. Rev. 
X, vol. 1, p. 021021, 2011. 

[52] H.-T. Wai, A. Scaglione, and A. Leshem, “Active sensing of social net¬ 
works,” to appear in IEEE Trans, on Signal and Information Processing 
over Networks, available online http://arxiv.org/abs/1601.05834v 1. 

[53] S. Boyd and L. Vandeberghe, Convex optimization. New York: 
Cambridge Univ. Press, 2004. 

[54] A. Shiriaev, Probability. Springer, 1996. 

[55] N. Friedkin, A. Proskumikov, R. Tempo, and S. Parsegov, “Network 
science on belief system dynamics under logic constraints,” under review 
in Science. 


Lemma 5: [44, Ch.XIII, §7]. A row-stochastic square 

matrix A is regular if and only if det(A/ — A) y 0 whenever 
X y 1 and |A| = 1; in other words, all eigenvalues of A except 
for 1 lie strictly inside the unit circle. A regular matrix is fully 
regular if and only if 1 is a simple eigenvalue, i.e. is the only 
eigenvector at 1 up to rescaling: Az = z ^ z = cld, c G R. 

In the case of irreducible [44] matrix A regularity and full 
regularity are both equivalent to the property called primitivity, 
i.e. strict positivity of the matrix A'" for some m > 0 
which implies that all states of the irreducible Markov chain, 
generated by A, are aperiodic [44]. Lemma 5 also gives a 
geometric interpretation of the matrix A*. Let the spectrum 
of A be Ai = l,A 2 ,...,Ad, where |Aj| < 1 as j > 1. 
Then can be decomposed into a direct sum of invariant 

d 

root subspaces R'^ = 0 Lj, corresponding to the eigenvalues 

i=i 

Xj. Moreover, the algebraic and geometric multiplicities of 
Ai = 1 always coincide [44, Ch.XIII,§6], so Li consists of 
eigenvectors. Therefore, the restrictions Aj = A\lj of A onto 
Lj are Schur stable for j > 1, whereas Ai is the identity 
operator. Considering a decomposition of an arbitrary vector 
V = y^j where Vj G Lj, one has A^ui = vi and A'^vj — >■ 0 
as fc —c» for any j > 1. Therefore, the operator A^ ■. v ^ vi 
is simply the projector onto the subspace Li. 

As a consequence, we now can easily obtain the equality 
(12). Indeed, taking a decomposition v = vi + ... + Vd, one 
easily notices that (/ — aA)~^vi = (1 — and (/ — 

aA)~^Vi —>■(/ — Ai)~^Vi as a ^ 1 for any i > 1. Hence 
lim (/ — aA)“^(l — a)v = vi = A*u, which proves (12). 


Appendix 

Stochastic regular matrices 


We state a spectral criterion for regularity. 
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